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SOLITONS ON A FINITE-GAP BACKGROUND 
IN BULLOUGH-DODD-JIBER-SHABAT MODEL 

I. YU. CHERDANTZEV and R. A. SHARIPOV 


Abstract. The determinant formula for TV-soliton solutions of the Bullough-Dodd- 
Jiber-Shabat equation on a finite-gap background is obtained. Nonsingularity condi¬ 
tions for them and their asymptotics are investigated. 


1. Introduction 

Found in the early works on an inverse scattering method [1, 2] the Bullough- 
Dodd-Jiber-Shabat equation 1 

u xt = e u -e~' 2u (1.1) 

makes a good alternative to the well-known Sine-Gordon equation as one more 
model of the two-dimensional integrable relativistic field theory with the self-action. 
In [3] the authors had suggested the construction of the finite-gap solutions for this 
equation. These are periodic and/or almost periodic function in x,t associated 
with special classes of double sheet ramified coverings of Riemann surfaces and 
represented by explicit formulae in terms of Pryrn theta-functions of such coverings. 
In this paper we construct Hirota type (see [4]) determinant formula for A-soliton 
solutions of (1.1) and use it for the investigation of the space-time asymptotics of 
them in space-time variables where (1.1) has a form 

u t t - u xx = e u - e~ 2u . (1.2) 

It is known that Eqs. (1.1) and (1.2) have no fast-decreasing soliton solutions being 
the same time real-valued and nonsingular. The presence of finite-gap background, 
as shown below, makes it possible for nonsingular solitons to decrease rapidly to that 
background. This has a transparent explanation in connection with the topology 
of appropriate Riemann surface. 

2. Main Construction 

First we touch briefly on the construction of finite-gap solutions of (1.1) being 
the background in the following considerations (for more details, see [3]). Equation 

1 S. P. Tsarev discovered that this equation was first found by Tzitzeica in [12]. Now it is called 
Tzitzeica equation. See also [13] for more details. 


Typeset by AmS-T^K. 



2 


I. YU. CHERDANTZEV AND R. A. SHARIPOV 


(1.1) has ’’zero curvature” representation [d x — L, dt — A] = 0 found in [5] with the 
matrix Lax operators of the form 


-U x 

0 

A 


0 

e~ 2u 

0 

1 

Ux 

0 

, 71 = 
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e u 
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1 

0 


A -1 e u 

0 
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According to the general scheme of finite-gap integration (see review [6]) one should 
find the Baker-Achiezer vector function xp = e(x, t, P) solving the spectral equations 

ipx = Lip, ipt = A%p (2.2) 

and depending on a point P of some compact Riemann surface T. As soon as 
matrices A and L are not of general form, the so-called “reduction problem” arises. 
In [3] we managed to solve it in our particular case in a way similar to that of [7] and 
[8] where the reduction of the general Schrodinger operator in magnetic field to the 
purely potential Schrodinger operator was found. Let T be the compact Riemann 
surface of the even genus g with the meromorphic function \{P) having the only 
pole Poo of the third order and the only zero P$ of the third order. In addition, let 
us suppose that T admits a holomorphic involution a with the property A(tr P) = 
—A(P) and the antiholomorphic involution r with property A (rP) = A(P). For this 
case, we can choose the local parameters k~ 1 (P) and g -1 (P) in the neighborhoods 
of marked points Poo and P$ defined by the following conditions k 3 = A, q~ 3 = A, 
k{rP) = k(P), q(rP ) = q(P). Baker-Achiezer functions satisfying (2.1) then are 
defined by fixing their essential singularities at Poo and Pq 


ipi 


k- 1 


ipi 


qe “ 

tp2 

~ 

k~ 2 

■ e kx , 

1p2 

~ 

q 2 e u 

V’3 


k~ 3 


^3 


q 3 


and by fixing the divisor D of their poles. Divisor D of degree g should satisfy the 
following limitations 


D + aD - P 0 - Poo ~ C, tD = D , (2.4) 

where C is a divisor of canonical class on F. The above limitations imposed on the 
choice of F and D make sure that the matrices L and A are of the form (2.1) and 
enable us to find the explicit formula for appropriate finite-gap solution u = v(x, t) 
of Eq. (1.1) in terms of the Prym theta-function qi(z) of the covering r F —> T/a: 

e v(x,t) =C ( F ) -2d xt lnn(Ux + Vt + z D ). (2.5) 

Vector zd is determined by the divisor D , constant c(r) depends only on Riemann 
surface T. Vectors U and V are determined by the normalized Abelian differentials 
of the second kind floe and Do, the main parts of them at Poo and Pq are of the 
form Dqo = dk + ... and Do = dq + ... respectively. 

Let r be the Af-curve with respect to anti-involution r (see [9]). In this case 
condition (2.4) is compatible with the choice of D such that each invariant cycle 
of r contains only one point of D except the cycle containing Poo and Pq where 
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there is no point of divisor D. It is the very choice which provides the real and 
nonsingular finite-gap solutions of Eq. (1.1) in (2.5). 

In order to construct the TV-soliton solutions of (1.1) let us fix an extra set of 
parameters consisting of 

(a) a set of numbers Ai, ... , Xn such that Af ^ Xj for i ^ j; 

(b) a set of nonzero constants C\ , ... , Cn- 

For general A,;, each equation A(P) = A,;, has exactly three solutions. Let us denote 
two of them as A,;, and <tA*, and then define the new fV-soliton Baker-Achiezer 
function *3 /(x,t,P) by the following requirements: 

(A) ^(x,t,P) is analytic everywhere on F except at P^, Pq and at points of the 
divisor D + Ai + ... + A n + A^ + ... + A']y; 

(B) it has essential singularities at P^ and P$ of the same form (2.3) as a purely 
finite-gap Baker-Achiezer function; 

(C) divisor V = D + Ai + ... + A n + A) + ... + is its divisor of poles and 

Res (3 vEq(P) A 2 (P) w(P)) = Cj ■ ^(aA*), 

Res (3 **(P) A 2 (P) w(P)) = -Cj ■ ^ t (aA J ). 

Here u>(P) = x>o,oo(P ) is an Abelian differential of the third kind with the divisor of 
zeros and poles D + crD — P 0 — Poo and with the residue +1 at Pq and the opposite 
residue at P^. The above conditions (A)-(C) fix up the function Baker-Achiezer 
ty(x,t,P) uniquely. From (2.6) we have 


N 

J2 A Res (^i(P)T 2 (ctP)A(P)w(P)) = 0, 


N 

J2 Res (^ 3 (P) d' 3 (crP) A 2 (P) w(P)) = 0. 

• i AjAj, 


This is the duality condition for U>(P) and \1/(<rP) which provides the equations 
(2.1) for fV-soliton Baker-Achiezer functions. 

In order to find the explicit formula for 7V-soliton solution of Eq. (1 1) we shall 
use the technique developed in [10] and [11], For each pair of vector functions ip 
and <p we define the pairing (ip | <p) as follows: 


(ip\<p) = -ipi<p2 +ip2(pi + X(P)ip 3 <p 3 . (2.7) 


The following properties of the pairing (2.7) are checked up immediately: 

d x (e(P) | e(aQ)) = (A(P) - A (Q)) e 2 (P) e 3 (P), 

a«(e(P)| eM »=e.(l-W) 


e 3 (P)ei(P). 


( 2 . 8 ) 
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Because of (2.8) the value of (e(P) | e(crP)) does not depend on x and t. One can 
calculate this value explicitly: 

(e(P)\e(aP)} = ^-. (2.9) 

3 u) 

Consider the following function depending on x, t and on the pair of points P and 
Q on Riemann surface T: 


f 1(x,t, P, Q) 


(e(x,t,P) | e(x,taQ)) 
A (P) ^ HQ) 


( 2 . 10 ) 


When Q is fixed function f2 has the pole divisor D + Q and essential singularities 
at P 0 and P QQ of the form 


fi(P,Q) 


q 2 e qt ei(aQ) 
A(Q) 


as oo, 


I2(P, Q)~ — k 4 e kx e 2 (aQ) as k —> oo. 


Moreover from (2.10) we derive 

Res (3 Cl(x,t, P, Q) A 2 (P) w(P)) = 1. 

P—Q 


( 2 . 11 ) 


( 2 . 12 ) 


Properties of 12 (x,t,P,Q) just mentioned give us the opportunity to use it for 
finding one of the components of 7V-soliton Baker-Achiezer function II /3 (x,t,P). 
We ’P 3 (x,t,P) construct via the following ansatz 


^3 (x,t,P) 


N 

= e 3 (x,t,P) + ^2 &{x,t,P,Aj) ■ aj(x,t) + 
l=i 

N 

+ ^2n(x,t,P,A*) ■ a*(x,t) 
l=i 


(2.13) 


with the parameters ay, a* yet undefined. We shall define them with the use of 
(2.6) by applying technique from [11]. 7V-soliton solution u(x,t ) of Eq. (1.1) then 
is found on a base of equivalences 


d t e 3 (P) - e v k ~ 4 e kx , d t ^ 3 (P) - e“ k ~ 4 e kx 

as P-> Poo. Taking (2.11) into account, one gets 

( N N \ 

e 2 (aAj) -aj+Y ^{crA*) ■ a* J . (2.14) 

l=i l=i / 

In general, to find aj and a* from (2.6) we should solve some system of linear 
differential equations, the special form of sums on the right part of (2.14) however 
makes it possible to eliminate this step and leads us to the formula 

e u = e v — d xt lndet(l — 12 C). (2.15) 
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Here C = diag(C 1 ,... , Cjv, — Ci,... ,—Cn) is a diagonal matrix matrix built 
of constants Cj from (b) (see above). Matrix O is defined by fl(x,t,P,Aj) and 
Cl(x, t, P , A*) with P = aClj and P = erA*. Therefore its components could be eval¬ 
uated explicitly via Riemann theta functions. This matrix is composed of blocks 



Bij = Cl(x,t,aAi,Aj) 

Cij = Q,(x,t,ak*,Aj) 

Dij = fl{x, t, crA*, A* 


Formula (2.15) with the matrix (2.16) is a direct analog of Hirota’s determinant 
formula from [4]. 

3. Problem of Nonsingularity 

In order to get real and nonsingular N-soliton solutions in (2.15) one should set 
more restrictions when choosing parameters in (a) and (b). Let the Riemann surface 


Fig. 1. The quarter of the Riemann surface 
r unfolded onto the plane, r acts as a reflec¬ 
tion downward, a acts as a rotation around 
the point P^. 

r and divisor D be chosen so that they correspond to the real and nonsingular 
finite-gap solution v(x,t) of Eq. (1.1). In this case surface T can be cut and then 
unwrapped onto a plane as shown on Fig. 1. Points where A(P) is purely real are 
shown with solid lines while purely imaginary values of A(P) are shown with dashed 
lines. A set of A(P) values on canonic a-cycles (i. e. invariant cycles of r except 
cycle ag containing Pg and Poo) consists of g intervals. When A j are in that interval 
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the canonic choice Aj, A* pairs preserving A(A j) = \(crA*) > 0 is possible. Each 
pair Aj, A* together with the point Pi, of divisor D define the orientation ij on a 
cycle a,i, in the direction from Aj to crAj through Pi. Thus we can prescribe the 
sign to Abelian differential w(P) with respect to orientation ij just introduced. Let 
us choose constants Cj, from (b) being real with the signs defined by 

sig n (Cj) = sign(o;(P)|jj) I . (3.1) 

\p=Aj 

It is easy to check that this rule gives an opposite sign on an opposite cycle era,; 
being in accordance with (2.6). 

Lemma 3 . 1 . When sign rule (3.1) holds, the number of zeros of the N-soliton 
Baker-Achiezer function ’F(P) on each cycle ai is equal to the number of its poles 
on that cycle. 

Since \E'(P) is real on a.j, (2.6) and (3.1) and simple geometrical considerations 
yield us the statement of lemma. This means that all zeros of 'F(P) are on a-cycles, 
hence they are separated from P 0 and P^, i.e. they never coincide nor come close 
to those points. The last fact in turn leads us to the following theorem. 

Theorem 3 . 1 . The above canonical choice of Aj, A* and the sign rule (3.1) for 
real constants Cj are enough to provide nonsingular N-soliton solution of ( 1.1) in 
(2.15). 

Such facts are well-known in the theory of finite-gap integration (see [9]). In the 
case of fast-decreasing solitons T has a topology of sphere with the only cycle ao 
containing Po and P^, and as a result we can find no place on T to set Aj and A* 
as required by the theorem above. 

Formula (2.15) together with the following estimates for purely finite-gap Baker- 
Achiezer function e(x,t,P) 

\e(x,t,P)\ < const • exp(/c oc (P) x + kq(P) t) (3.2) 

via the Abelian integrals 


Koo (P) = Re 



kq{P) = Re 



fio 


give us the opportunity to find asymptotics for TV-soliton solutions of Eq. (1.1) as 
x.t —> ±oo. Because of (3.2) we may divide entire rows and columns of matrix Q 
in (2.15) by exponential factors of growth 

exp(2 (ko^ctA*) - Koo(aAj)) x + 2(k 0 (ctA^) - n 0 (crAj))t). 

Vanishing of (^(ctA*) — K oc (aAj)) x + 2(k 0 (ctA*) — Ko(<rAj)) t determines the free 
soliton trajectories. Their velocities are 

Koo(cA *) - Koo(aAj) 

Vj Ko(aAj) — Ko(aAj) 


(3.3) 
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For Eq. (1.2) in space-time coordinates soliton velocities Vj, are bound with Vj as 


Vj = 


e £ Vj + e E 

e -s — e £ 


(3.4) 


Here £ is a scalar parameter depending on the choice of particular system of 
Lorentzian coordinates. The matter that Vj are less than light velocity c = 1 
here is a consequence of (3.3) and (3.4) and the following lemma. 

Lemma 3.2. When points A j, a A* are placed as shown on Fig. 1, the inequalities 
K 00 (aA*) — Kao(aAj),0 and ko(<jA*) — no(aAj) < 0 hold. 

The observers being at rest in almost all Lorentzian systems of coordinates will 
see the similar pictures of soliton interaction: before and after the interaction soli- 
tons are free. They are separated by finite-gap background potentials v(x,t,Zj) 
with the common spectrum T but different phases Zj , peculiar to each interval be¬ 
tween solitons. The background phase shift obtained while crossing the path of j-th 
soliton is A Zj = U (A j + A* — crAj + a A j). Here A Zj is a vector of Prym variety of 
the covering T —> V/a and U is a map from subset of Div(T) into Prym(T) built in 
paper [3]. 
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